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ABSTRACT 

A photon emitted near a compact object at an angle a with respect to the radial direction escapes to infinity 
at a different angle ip > a. This bending of light is caused by a strong gravitational field. We show that, in a 
Schwarzschild metric, the effect is described by 1 — cos a = (1 — cos^)(l — Tg/R) where R/rg is the emission 
radius in Schwarzschild units. The formula is approximate and it applies at i? > 2rg only, however at these radii 
it has amazing accuracy, fully sufficient in many applications. As one application we develop a new formulation 
for the light bending effects in pulsars. It reveals the simple character of these effects and gives their quantitative 
description with practically no losses of accuracy (for the typical radius of a neutron star R = 3rg the error is 
1%). The visible fraction of a star surface is shown to be Sv/AttR'^ = [2(1 — Vg/R)]^^ which is 3/4 for R ~ 3rg. 
The instantaneous flux of a pulsar comes from one or two antipodal polar caps that rotate in the visible zone. The 
pulse produced by one blackbody cap is found to be sinusoidal (light bending impacts the pulse amplitude but not 
its shape). When both caps are visible, the pulse shows a plateau: the variable parts of the antipodal emissions 
precisely cancel each other The pulsed fraction of blackbody emission with antipodal symmetry has an upper 
limit Amax — {R — 2fg)/(i? + 2rg). Pulsars with A > Amax must be asymmetric. 

Subject headings: gravitation — radiation mechanisms: general — relativity — stars: neutron — (stars:) 
pulsars: general — X-rays: binaries 
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1. INTRODUCTION 

Light bending by gravity is a classical effect of general rela- 
tivity (e.g. Misner, Thorne, & Wheeler 1973). In a spherically 
symmetric gravitational field, the exact bending angle is given 
by an elliptic integral (see Appendix). The effect plays impor- 
tant role in astrophysics; e.g. weak bending gives rise to the 
gravitational lensing phenomenon. Strong bending occurs near 
compact objects (black holes and neutron stars) and crucially 
affects their emission. 

In this Letter we give a simple formula for the bending an- 
gle that replaces the elliptic integral with high accuracy at radii 
R > 2rg. Here Vg = 2GM/c^ is the Schwarzschild radius of 
the gravitating spherical object of mass M. We apply this for- 
mula to neutron stars and develop a new simple formalism for 
the light bending effects in pulsars. We illustrate with pulsars 
that have two point-like antipodal spots emitting thermal radia- 
tion, and compare the results with the exact theory (Pechenick, 
Ftaclas, & Cohen 1983). The advantage of our formalism is 
that it gives a clear understanding of the bending effects and 
describes them analytically with almost no losses of accuracy 
as long as stars with R > 2rg are considered. The exact theory 
well describes stars with any R, including R < 2rg where the 
bending angle (3 can exceed 90°, however it requires a compli- 
cated numerical treatment. Standard models of neutron stars 
predict R > 2rg with a typical R w 3rg (e.g. Shapiro & 
Teukolsky 1983). Our formalism then works with ~ 1% ac- 
curacy and should be sufficient in many practical problems. It 
should also be good for accretion disks around Schwarzschild 
black holes: the disk inner edge is at 3rg and the bulk of emitted 
radiation does not enter the region R < 2rg. 



2. THE COSINE RELATION 

Let an emitter be located at point at a radius R (Fig. 1). 
A photon emitted along the radial direction escapes radially. 
A photon emitted at some angle a with respect to radius (a 
is measured by the local observer at E) escapes along a bent 
trajectory. We will describe the trajectory in polar coordinates 
(r, ip) where = for the escape direction. The standard way 
of computing a{il)) is given in Appendix. Instead one can use 
the approximate relation 



cos a 



(1-cos^) (l-^). 



(1) 



The accuracy of relation (|l]) is remarkably high. It gives the 
bending angle /? = t/j— a with a small fractional error e — 6(3/ fi 
shown in Figure 2. For example for R = ivg the maximum er- 
ror is e,„ax = 3% (at a = 90°) and e < 1% for a < 75°. High 
accuracy is maintained for R > 2rg. For R < 2rg the approx- 
imation is not applicable: then cos a ~ would correspond to 
cos?/) < —1. We therefore limit our consideration to radiation 
escaping from R > 2rg where (3 < 90°. 

The high accuracy of equation (|l]) is a striking property of the 
Schwarzschild spacetime. This functional shape cannot be un- 
derstood as a linear expansion in rg/R (that expansion is given 



in Appendix and its accuracy is much worse). Equation (|1|) can 
be understood as follows. Let us consider a; = 1 — cos a as 
a small parameter and expand y = 1 — cos in terms of . 
After some algebra we find 
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where u = TgjR. The linear term is just equation ^\ it de- 
scribes the standard solid angle transformation d£lQ — 2Trdx 
dil — 2TTdy along the radial ray. The interesting property of 
equation (0) is that the x^-term vanishes and the higher order 
terms have small coefficients. It explains the extremely high ac- 
curacy of equation (|lj) at small u and/or x (Fig. 1). Only when 
both a; ^ 1 and u 0.5 there is a sizable correction to the 
linear term and here e peaks sharply. The linear expansion in u 
(also derived in Appendix) reads y — x{l + u) + 0{v?). 

Though the cosine relation (|l]) is sufficient for our subsequent 
discussion, note that it also gives a simple description for the 
shape of a photon trajectory with a given impact parameter h 
(Fig. 1). Equation ([l|) can be written at any point (r, V') along 
the trajectory if a measures the angle between the photon ve- 
locity and the local radial direction. Combining ([l|) with dlT 
we find = 6^[sin^ ?/) + (1 — cos V-')^?'g/?' 



and 



r1{l — cosil^Y 6^ 

— 1 

4(1 -h cos ■0)^ sin^i/) 



1/2 



rg(l — COS V") 
2(1 + cos?/') ' 



(3) 



Trajectories (||) are very close to the exact ones. The fractional 



error of r at a given is (5r/r < 4 x 10 



as long as r > 2rg . 



3. THE VISIBLE SURFACE OF A STAR AND THE OBSERVED FLUX 

Consider a spherical compact star viewed by an observer at 
distance D R. Choose spherical coordinates {R, i/j, (f) with 
the polar axis = directed toward the observer. What is the 
observed flux dF from a surface element dS of the star? First 
of all one should decide whether dS is visible. In a flat space- 
time the visibility condition would be simply /i = cos t/; > 
where is the polar angle of dS. The star gravity creates cur- 
vature that allows the observer to see regions with negative jj. 
down to a critical jiy that defines the dark side of the star The 
observed light from the circle ji = /i„ is emitted tangentially to 
the star {a = 90°). Exact calculations of a{'4)) via the elliptic 
integral yield tabulated /i^,, e.g. /i„ = -0.886, -0.633, -0.484 
for R/rg = 2, 2.5, 3. From equation ([l|) we find immediately 
the approximate ji^ — —rg/{R — rg). The visible fraction of 
the surface is Sy/AnR^ = (1 - Mi,)/2 = [2(1 - Tg/R)]-^, e.g. 
we get 0.75 for R — 3rg while the exact value is 0.742. 

Appendix gives the standard formula for the observed flux 
from a visible dS (eq. [ll]). It requires the knowledge of a{ip) 
and the exact theory leads to a complicated numerical treatment 
of the problem (Pechenick et al. 1983). Instead we use the ap- 
proximate equation (pi); then equation (pll) reads 



^-^=(1^"B") -^o(a)cosa^, 



R 



(4) 



where cos a ~ — Vg/R) + rg/R- Equation (Q) is surpris- 
ingly simple. Apart from the redshift factor (1 — Vg/R)'^ the 
element dS emits as if the problem were Newtonian with ip re- 
placed by a. 




Fig. 1 . — Example photon trajectory. The emission point E isai R = 2.5rg 



and the emission angle is cos a 



0.1; 



eter b = 3.21rg. In this case cos?/) 



the corresponding impact param- 
0.482L and the bending angle 



/3 = i/i — o 35°. The approximate trajectory (eq. ^ is also plotted here and 
it is not distinguishable from the exact trajectory. 
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Fig. 2. — Accuracy of the cosine relation Alh. This relation gives the bendig 
angle (5 = xj) — a with an error e = (5/3/p which depends on the emission 
radius R/vg and the emission angle a. Here the contours of e = const (with 
logarithmic step 0.5) are plotted on the TgjR — a plane. 



When the local intensity of radiation /o(<^) is known for 
all dS, equation (^ can be integrated over the visible surface 
fi > fiv get the total observed flux. For example, suppose a 
star has a homogeneous hot spot of arbitrary size and shape that 
emits thermal radiation [Io{a) = const]. Denote the spot area 
by s = j dS and the corresponding projected area by s± = 
J fidS (if only a part of the spot is in the visible zone fi > fj,^ 
the integrals are taken over this part). Then the observed flux is 
F = {Io/D^){l - rg/R)^[s±{l - rg/R) + srg/R]. Note that 
< when the spot is in the region /i^, < < 0. 
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4. PULSARS 

Consider a neutron star with two antipodal hot spots associ- 
ated with polar caps of the star magnetic field. Assume here 
that the spot size is small compared to the star radius R. The 
star is spinning and its magnetic axis is inclined to the spin axis 
by an angle 9 < 90°. Therefore the spots periodically change 
their position and inclination with respect to a distant observer 
and the observed radiation pulsates. 

The spot whose circle of rotation is closer to the observer will 
be called "primary" and the other symmetric spot "antipodal". 
Denote the unit vector pointing toward the observer by d and 
the angle between d and the spin axis by i. To compute the ra- 
diation fluxes from the primary and antipodal spots one needs 
to know their inclinations: /i = n • d and ft = n ■ d, where n 
and n — — n are the spot normals. When the star executes its 
rotation, /i varies periodically between /imin — cos(« + 9) and 

Mmax = COs(i - 6), 



IJ,{t) = sin 6 sin i cos fit + cos ( 



cos z. 



(5) 



Here is the angular velocity of the pulsar; t — is chosen 
when n = ^max- 

Consider the blackbody component of the pulsar emission; it 
has isotropic intensity Io{a) — const. Using we find the 
observed fluxes F and F from the primary and antipodal spots 
(when they are visible) 



F 



/i 1 



R' 



F 



-Ml 



R' 



(6) 



Here Fi = (1 — rg/R)'^lQs/D'^, s is the spot area, and point- 
like-spot approximation is adopted. For the antipodal spot we 
substituted p, — —fi. The observer sees the primary spot when 
/i > /i„ and the antipodal spot when fl > /i„ (§ 3). De- 



note K 



\^v\- Both spots are seen when 



-K < fl < K and then the observed flux is 



-Fobs — F 



F = 2^F, 



const ! 



(7) 



Hence the blackbody pulse must display a plateau whenever 
both spots are in the visible zone. We emphasize the high accu- 
racy of equations (||) and (Q), though the exact theory does not 
give the ideal plateau. The fractional error 6 in Fobs is maximal 
when n = e.g. Jmax ~ 1.3% fori? = 3rg. For R = 2rg our 
formalism predicts that the whole surface is visible and Fobs is 
given by equation (0) at any 9,i, fi; then (Smax ~ 6%. 

At any time Fobs takes one of three values: F, F, or F + F. 
A pulsar shows the plateau if cos(i + 9) < k. A pulsar 
with arbitrary i and 9 belongs to one of four classes described 
in Figures 3 and 4. It is easy to compute the pulsed fraction 
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A reaches its maximum Amax if the pulsar has fin 
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< K. This maximum is ^max — {R~ '^fg)l{R + 2rg). 
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Fig. 3. — Location of classes I-IV on the i — 8 plane. The angle f is 



defined so that cos ^ 



g), e.g. S, = 60° for R 



= 3rg. The 
cos(i + 0), 

> K. 



classes are determined by the relative positions of 
Mmax = cos(j - 6), K = rg/{R - Vg), and 
The antipodal spot is never seen in such a pulsar (and the piimary spot is visi- 
ble all the time). The observed blackbody pulse has a perfect sinusoidal shape. 
Class II: —K < ^miii < K < ^max- The primary spot is seen all the time and 
the antipodal spot also appears for some time. The pulse has a sinusoidal pro- 
file F{t) interrupted by the plateau m. Class III: 

Mmin < — z^- The primary 
spot is not visible for a fraction of period (and then only the antipodal spot is 
seen). Suchj pulsar shows a primary sinusoidal profile F{t) interrupted by 
the plateau ( 71) and the plateau is interrupted by a weaker sinusoidal sub-pulse 
F{t) from the antipodal spot. The sub-pulse occurs when only the antipodal 
spot is in the visible zone. Class IV: ~k < fimim Mmax < Both spots are 
seen at any time. Then the observed blackbody flux is constant (eq. M). Only 
anisotropic emission with Io{c() const can contribute to pulsations. 
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Fig. 4. — Blackbody pulse profiles for classes I-IV defined in Fig. 3. 
R = 3rg is chosen and the flux is normalized to F\ (the maximum possi- 
ble flux that is observed when the primary spot is viewed face on). Class I is 
exemplified with Q = 20° and i = 30°, class II with 9 = 30° and i = 60°, 
class III with 9 = 60° and i = 80°, and class IV with 9 = 20° and i = 80°. 
The approximate model developed in this paper is shown by solid curves and 
the exact model is shown by dotted curves. The plateau has the same level in 
all cases: Fp/Fi = 2rg/R. = 2/3. 
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5. CONCLUSIONS 

A new formulation for the light bending effect is given in 
this Letter It makes use of an approximate simple formula that 
describes photon trajectories with remarkable accuracy (eq. 1). 

This formulation is applied to pulsars where light bending 
is known to be a crucial effect. The developed analytical the- 
ory replaces the previously used complicated numerical treat- 
ment with high accuracy, sufficient for practical purposes. For 
example, it immediately gives the observed pulse profile for a 
rotating neutron star with two antipodal hot spots (eqs. |] and 
^). These analytical expressions highlight the effects of light 
bending in pulsars: (i) If the antipodal spot is never seen, the 
blackbody pulse is perfectly sinusoidal. The bending affects 
the pulsed fraction A only (eq. (ii) The visible surface of the 
star is significantly enhanced, S^/AttR^ — [2(1 — rg/R)]~^, 
and both spots may be seen for a fraction of period or all 
the time. Then the blackbody pulse shows a plateau Fobs — 
Fp — {2rg/R)Fi where Fi is the flux from one spot viewed 
face on. (iii) The observed flux cannot exceed Fi and it can- 
not be smaller than Fp. The latter is a purely relativistic ef- 
fect: Fp = in the Newtonian limit rg/R 0. The max- 
imum pulsed fraction of a pulsar with antipodal symmetry is 
Anax = (i^i - Fp)/{Fi + Fp) = {R - 2rg)/{R + 2rg). This 
is a severe constraint, e.g. A^^ax = 1/5 for R = 3rg instead of 
Newtonian A^^ax = 1- 

Within a few percent the results agree with the previous nu- 
merical works. The new formalism, however, is much simpler 
in use. It gives a clear description of the pulse profiles for given 
i, 6, and R/rg and makes it easier to compare the theory with 
observations and constrain pulsar parameters. We focused in 
this Letter on pulsars with small hot spots with opening an- 
gles 69 < 10°. The extension to big spots is straightforward 
(§ 3). Note here that the plateau is predicted when both spots 
are fully visible and the upper limit ^max will remain valid for 
blackbody emission with antipodal symmetry. 

Caps of real pulsars (both radio and X-ray) do not emit as 
a blackbody and their Iq (a) need not be isotropic. In partic- 
ular: (1) X-ray pulsars show Comptonized power-law spec- 
tra, sometimes along with a clear blackbody component (as in 
SAX J1808.4-3658). Light bending affects the two components 
differently: the blackbody pulse from one spot remains sinu- 
soidal while the Comptonized pulse is distorted. This distortion 
is described by equation (Q) and it depends on the anisotropic 
intensity Io{a) of the Comptonized emission. (2) For pul- 
sars of high spin (with ms periods) the Doppler effect changes 
Io{a). It is easily added to the model by the corresponding 
transformation of Iq from the spot rest frame to the local non- 
rotating frame. (3) A strong magnetic field B ~ 10^'^ G leads 



to anisotropic emissivity with a preferential direction along the 
field. In any case, calculations of Io{a) separate from the light 
bending effects. Radiation that propagates freely through the 
pulsar magnetosphere is described by equations ([l|) and (Q) for 
any /o(a). 

APPENDIX: LIGHT PROPAGATION IN A 
SCHWARZSCHILD METRIC 

For a given photon trajectory let us choose Schwarzschild co- 
ordinates x'' — {t, r, 9, tp) so that the trajectory is in the plane 
9 — 90°. Let u'^ = dx'^/dX be the 4-velocity of the photon 
(A is an affine parameter, A = r at r ^ oo). The trajectory 
has integrals ut = 1 and — b. From u^Ui — one finds 
{u^y = 1 — (5^/r^)(l — Tg/r). Let us choose ^ = for the 
escape direction. Then at radius R the photon has 

(9) 

The emission angle a (between the emitted photon and the local 
radial du^ection) is tan a = ^Ju^u^/ y^uTu^ which yields 



sin a = 




Combining equations (^ and ( |1C| ) one can compute numeri- 
cally the relation between ij} and a for a given R. At small u = 
Tgl R one can expand equation ^ in u and keep the liner term 

only. Using the equality /^'""(sina — z'^){l — z'^)^'^/'^dz = 
2(1 — cos a) one gets tp — a + u{l — cos a)/ sin a + O(u^). 

Consider a star of radius R with a surface radiation inten- 
sity /o(a) which can vary over the surface. The observer lo- 
cation is given by a radius-vector D (I? ^ R). A (visible) 
surface element dS — R^dfi dcf) is observed at impact parame- 
ters (6, b + db) and subtends a solid angle dfl = bdb dcj)/ D"^ on 
the observer sky. Here /i — cos ?/; and is an azimuthal angle 
corresponding to rotation around D; b depends on fi only, not 
on (f). The observed radiation flux from dS is dF = Idfl where 
/ = {l — rg/R)^Io [/and/o are related through //t^* = lo/'^o 
and the frequency redshift I'/i'o = (1 ~ ^^e Misner 

et al. 1973]. Using equation ( p^ ) one gets 

, ^ lb .db ,dS ( r„\ _^ , . d cos a dS 
'^'^=R^^dil^D^ = V- r) 

(11) 
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